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Abstract. Particle-in-Cell simulations of magnetic reconnection with an 
H + current sheet and a mixed background plasma of H + and + ions are 
completed using physical mass ratios. Four main results are shown. First, 
the + presence slightly decreases the reconnection rate and the magnetic 
reconnection evolution depends mainly on the lighter H + ion species in the 
presented simulations. Second, the Hall magnetic field is characterized by a 
two-scale structure in presence of + ions: it reaches sharp peak values in 
a small area in proximity of the neutral line, and then decreases slowly over 
a large region. Third, the two background species initially separate in the 
outflow region because H + and + ions are accelerated by different mech- 
anisms occurring on different time scales and with different strengths. Fourth, 
the effect of a guide field on the + dynamics is studied: the + presence 
does not change the reconnected flux and all the characteristic features of 
guide field magnetic reconnection are still present. Moreover, the guide field 
introduces an + circulation pattern between separatrices that enhances high 
+ density areas and depletes low + density regions in proximity of the 
reconnection fronts. The importance and the validity of these results are fi- 
nally discussed. 
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1. Introduction 

The presence of heavy + ions with ionospheric origin has long been observed in Earth's 
magnetotail [Frank et al, 1977; Wilken et al, 1995; Zong et al, 1998; Mouikis et al, 2010]. 
+ population increases during storm times when + can even dominate H + in number 
density and pressure [Kistler et al, 2005] and there is evidence that + species is involved 
in magnetic reconnection events [Kistler et al, 2005; Wygant et al, 2005]. Besides finding 
the presence of + in the Earth's magnetotail during magnetic reconnection, the analysis 
of satellite data showed different behaviors of + and H + populations [Wygant et al, 
2005; Kistler et al, 2006]. However, pure data analysis can not yet fully explain the 
dynamics of H + and + because of the lack of resolution in space of satellite data. In 
support of observational analysis, theoretical and simulation studies can help first in 
verifying observations and second in unveiling the microphysics of magnetic reconnection 
in presence of mixed ion populations. In fact, simulations can be used as a microscope to 
study dynamics at micro-scale and to reveal the microphysics of magnetic reconnection. 

Many simulation approaches can be applied to study the effect + on the magnetosphere 
dynamics and on magnetic reconnection. Global simulations of the Earth's magnetosphere 
with + ionospheric outflow have been completed [Winglee et al, 2002; Brambles et al, 
2010a; Glocer et al, 2009a, b; Wiltberger et al, 2010; Brambles et al, 2010a, b]. The 
multicomponent fluid approach focused more specifically on the study of magnetic recon- 
nection with + ions [Shay and Swisdak, 2004]. Moreover, the simulation of test particles 
moving in a preassigned electromagnetic field, that is calculated previously with a MHD 
simulation, allowed to study + ion trajectories in the global Earth's magnetosphere [Birn 
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et al, 2004]. Hybrid simulations, where + and H + ions are modeled with kinetic parti- 
cles while electrons are treated as a fluid in a self-consistent field, have used also to study 
magnetic reconnection with + ions [Fujimoto and Nakamura, 1994]. Finally, a fully ki- 
netic approach, such as the Particle-in-Cell method, where electrons are represented by 
particles also, has been followed by Hesse and Birn [2004]; Karimabadi et al. [2010]. 

The present paper reports the first full kinetic simulations with the physical mass ratios 
among electrons, hydrogen and oxygen. Previous Particle-in-Cell simulations either did 
not use physical mass ratios for + , H + and electrons [Hesse and Birn, 2004; Karimabadi 
et al, 2010] or did not consider + ions in the magnetic reconnection dynamics [Ricci 
et al., 2002, 2004]. 

The aim of this paper is to study the effect of a background + population on the 
evolution of magnetic reconnection using realistic parameters. A current sheet composed 
of H + only, and a mixed background of + and H + , to mimic the lobe population of the 
magnetotail, are considered. The results are compared with identical simulations but with 
only H + plasma background. After an overall study of magnetic reconnection evolution, 
an in depth analysis of the flow pattern and on the different dynamics of background + 
and H + particles is performed to study the role of + in magnetic reconnection. 

The paper is organized as follows. The second section presents the simulation parame- 
ters and the initial set-up. Three main results are reported in the third section. First, the 
overall magnetic reconnection evolution and the electromagnetic field at representative 
times and the different acceleration mechanisms acting on H + and + ions are presented. 
Second, the density evolution of the two species is studied and the initial H + + sep- 
aration in the outflow region is shown. Third, the effect of a guide field on magnetic 
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reconnection with + and the H + dynamics is analyzed. The fourth section concludes the 
paper summarizing the results and discussing the validity of the presented simulations. 

2. Simulation approach and physical system 

Two and half dimensional (two spatial components and three components for the electro- 
magnetic fields and particles velocities) Particle-in-Cell simulations have been completed 
on the X-Y plane with X along Earth-Sun direction in the Earth's magnetosphere, and 
with Y pointing to the north-south direction. The X and Y coordinates correspond to 
the -X and Z coordinates respectively in the GSM system, while in our coordinate system 
Z is out-of-plane ignorable coordinate. An antiparallel magnetic field, vanishing in the 
mid-line of the simulation domain (y = L y /2), is initialized in the X direction as: 

B x = B tznh( y ~^ y/2 ) (1) 

where A is the current sheet thickness, and B Q is the lobe magnetic field. The number 
density has the following profile: 

n(y)=n cosh- 2 (^^) (2) 

where no is the peak particle density at the center of the system. The simulation box is 
AOdn (L x ) and 20d H (L y ) long in the X and Y directions, where dn is the H + skin depth 
based on the peak density n of the Harris sheet with A = Q.hd H . The particles forming 
the current sheet are electrons and only H + ions. They are initialized with a Maxwellian 
distribution with a velocity drift that satisfies the Harris equilibrium force balance. The 
electron temperature is (Te/fmgC 2 )) 1 / 2 = v t he/c = 0.045, where c is the speed of light in 
vacuum, m e and v t he are the electron mass and thermal velocity. The H + ions have five 
times the electron temperature, T H = 5T e . In addition, a tenuous background plasma 
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with H + and + densities ribH = nbo = 0.05no is added to mimic the lobe plasma. This 
background plasma is uniformly distributed in the simulation box; it is not drifting and 
has the same temperature of the current sheet plasma T H = T Q . The simulation time 
step is ojpnAt = 0.06, where oj p h = ^4:7in e 2 /mH is the H + plasma frequency, m H is 
the H + mass and e is the elementary charge. In this simulation set-up, the ratio u p h to 
Vl cH = B e/m H (H + gyro-frequency), that is equal to the ratio of c to Vah, the H + Alfven 
velocity calculated with B = B and n = n , results 275. This ratio is comparable to the 
typical value of 300, observed in the Earth's magnetotail [Lapenta et al, 2010]. The grid is 
composed of 1024 x 512 cells. In total 411,041,792 computational particles with physical 
mass ratios m#/m e = 1836 and mo/ ran = 16 are used in all the reported simulations. 
A localized perturbation around the X point is applied to initiate magnetic reconnection 
without a boundary driving avoiding the creation of spurious compressive waves in the 
simulation domain [Lapenta et al, 2010]: 

5A Z = A z0 cos(27r(:r - L x /2)/L A )cos(7i(y - L y /2)/L A )e^ L ^ 2+ ^- L y/ 2 ^ 2 , (3) 

with La = lOcr and a = 0.5dn- The boundary conditions are periodic in the X direction, 
while conducting and reflecting boundary conditions are imposed in the Y direction for 
the electromagnetic field and particles. A simulation with only H + ions in the background 
plasma has been completed to compare the results and determine the effect of the + 
species in the magnetic reconnection. The companion simulation has identical parameters 
but starts from a background population, that is composed only of H + with initial density 
nbH = 0.1n . Simulations are completed using the parallel Particle-in-Cell code iPICSD 
[Markidis et al, 2010], that uses the moment implicit Particle-in-Cell scheme and allows 
large time step and grid spacing still retaining the numerical stability. 
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3. Simulation results 

The reconnected flux is calculated as the difference of the magnetic flux between 
the X and O points and it is studied to characterize the speed of magnetic reconnection. 
Its evolution is shown for the two simulations with a mixed population of + and H + ions 
and with only the H + species in panel a) of Figure 1. A slow growth of the reconnected 
flux is visible until Vt cH t « 5, then the fast phase starts. The + presence slightly affects 
the evolution of magnetic reconnection: the reconnected flux is approximately the same 
for both simulations. The reconnection rate is calculated as the derivative in time of the 
reconnected flux computed every 1.1VL~^. It is normalized to the BqVah/c, and shown in 
panel b) of Figure 1. It increases until time fl c nt ~ 10, and reaches the peak value ap- 
proximately at time fl c nt ~ 12. When the two simulations with different ion background 
are compared, the reconnection rate peaks differ by 0.034 B Vah/c, approximately 10 % 
of the maximum reconnection rate. 

The pseudocolor plots of different components of the electric field at time fl C Ht = 10.9, 
when the reconnection rate is maximum, are presented in three panels of Figure 2. The 
reconnection electric field, E z in the chosen reference system, is shown in panel a). The 
peak value of the reconnection field is 8 x lO~ 6 (uj pH m H c)/e. The in-plane electric fields 
are localized along the reconnection fronts, and are characterized by higher values around 
2 — 4 x 10~ 5 (u P HiTiHc)/e, as visible in the pseudocolor plots in panels b) and c) of Figure 
2. 

A comparison of the quadrupolar out-of-plane component of magnetic field, signature 
of Hall reconnection, at time VL c nt = 16.4 is shown in Figure 3. The B z component for 
the simulation with and without the + background species are presented respectively 
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in panels a) and b). Clearly, the typical quadrupolar structure reaches the same peak 

values (approximately 0.0009oj p Hm H c/e = 0.25B ), but it is broader in the case of the 

simulation with mixed H + , + background populations. In addition, Figure 4 presents 

the line profile of the out-of-plane component of the magnetic field B z along x = 8dn to 

compare the peak values and the thickness of the quadrupolar structure. In proximity of 

the neutral line (y = lOdn), two small B z peaks reveal the electron scales. Moreover, the 

Hall magnetic field profile in presence of + (red line in Figure 4) reaches peak values in a 

relatively small region, starting from the neutral line and ending at 4dn from it, and then 

slowly decreases over an extended area reaching the boundary. Thus, the B z profile can 

be divided in two parts: the sharp peaks representing the "core" Hall magnetic field, and 

a slowly decreasing B z representing the "tail" Hall magnetic field. This two-scale nature 

of the Hall magnetic field reflects the presence of the two H + and + dissipation layers, 

and of light and heavy whistler scales [Shay and Swisdak, 2004]. 

The magnitude of the average magnetic force divided by the charge e (v s x B/c), where 

v s is the average velocity for species s, is plotted in Figure 5 to compare the effectiveness of 

the electric and magnetic fields particle acceleration on H + and + species. It is clear from 

panel a) of Figure 5 that the magnetic force acting on the H + background population is 

maximum along the separatrices and in a shell around the X point. The immediate region 

around the X point has vanishing magnetic field, and the H + ions in it are demagnetized. 

The maximum intensity value is 8.42 x lQ~ & ui pH mHc/ e. The average magnetic force acting 

on + particles is shown in panel b), and its maximum values is four times less than the 

H + one, 2.14 x 10~ 6 u pH m H c/e. 
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3.1. Acceleration mechanisms 

It is useful to study the flow patterns and velocities for the two ion background species to 
understand different acceleration mechanisms. Figure 6 shows the in-plane (X-Y plane) 
average velocity for the H + and + background populations at time Vt cH t = 18.5. At 
this time, H + background ions have high velocity mainly in the outflow region. The H + 
peak velocity reaches the 0.00461c value (1.27Vah)- On the other hand, the + species is 
characterized by a rather different velocity pattern as clear in panel b): its flow is localized 
in the inflow region, and it is normal to the reconnection separatrices. The peak velocity 
is 0.001405c (0.39Vah), approximately three time less than background H + velocity. In 
addition, a small part of the + population in the X point region participates in the 
outflow with an approximately 0.0005c velocity (Q.IAV AH )- 

A plot of the streamlines at time f2 c #t = 18.5 is shown in Figure 7. The streamlines 
are calculated as tangents to the average velocities in the X-Y plane at a given time. In 
panel a) the H + flow, represented by black lines, moves from the inflow region towards 
the X point and it is deflected in the outflow region. The thickness of the outflow is 
approximately 3d#. Instead, the + ions follow a rather different flow pattern, which 
is represented with red lines in panel b). A large part of the + ions moves from the 
outflow and is deflected perpendicularly to the separatrices in proximity of the reconnec- 
tion fronts. A small part of the + population is deflected by a 90° angle around the X 
point. The simulations show that the acceleration acting on + ions could be explained 
by the mechanism proposed in Refs. [Drake et al, 2009b, a]: + particles behave like 
pick-up ions and are accelerated by magnetic reconnection exhausts in the boundary layer 
between the upstream and downstream plasmas. 
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3.2. Species separation 

The different acceleration and dynamics of background ions have important macroscopic 
consequences. The H + background population is accelerated towards the X point by the 
E x B drift [Pritchett, 2001] and expelled in the outflow region. The H + species remains 
concentrated on the dipolarization fronts. Instead, the + component does not effectively 
experience a magnetic drift and responds mainly to the normal electric field in corre- 
spondence of separatrices. These different acceleration mechanisms determine enhanced 
density areas of H + , and + in different regions. Figure 8 shows a pseudocolor plot of 
the density, normalized to no, at time fl C Ht = 18.5 for the + and H + background pop- 
ulations. The initial value density for both species is ribu = ^f>o = 0.05n . A comparison 
between the two panels of Figure 8 reveals that peak density is higher in the case of H + , 
that reaches a value of approximately 0.2n , while the peak density for + population is 
0.15n . Moreover, the two background species concentrate in different parts of the sys- 
tem. This is clear from panel a) of Figure 9, a pseudocolor plot of the peak densities at 
0, cH t = 18.5. The peak density plots show the regions with density values between the 
maximum and 70% of the maximum. The background H + , + ions and electrons appear 
respectively as black, red and green areas. In particular, the H + and electron peak density 
is localized on the dipolarization fronts. Instead, the + population is slower and follows 
the H + species in a more extended region from the tail of the H + peak density to the X 
point. Panel b) of Figure 9 shows the line profile of the three species densities along the 
y = L y /2 line using the same color scheme of panel a). Figure 10 presents the H + and 
+ peak densities at three different times, showing the progressive separation of the two 
ion species in the outflow region at the initial stage of magnetic reconnection. 
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3.3. The effect of a guide field 

To mimic the configuration of magnetic reconnection in the magnetopause, it is useful 
to add an uniform out-of-plane guide field in the initial Harris sheet equilibrium [Ricci 
et al, 2004; Swisdak et al, 2005; Lapenta et al, 2010; Goldman et al, 2010; Lindstedt 
et al, 2010]. A guide field, perpendicular to the X-Y plane, B g = B is imposed in the 
initial magnetic configuration. Figure 11 shows a comparison between the reconnected 
flux between two simulations with + and with and without the guide field in red solid 
and dashed lines, and a simulation with only background H + and guide field in black 
line. Both simulations with guide field present a slow phase where the reconnected flux 
increases slowly, that is followed by a fast phase. The reconnected flux is almost identical 
for both simulations with and without + in presence of a guide field. By comparing the 
reconnected flux with and without the guide field, it is clear that a guide field B reduces 
the speed of magnetic reconnection to approximately 32% in the case of simulations with 
only H + ions. In fact, the reconnection rate peaks for simulations with and without 
guide field are respectively 0AB Vah/c and 0.27B Vah/c These results are in reasonable 
agreement with those reported by previous studies [Pritchett, 2001; Ricci et al, 2004]. 

Panels a) and b) of Figure 12 show a quiver plot of the average velocity for background 
H + and + populations at time VL cH t = 21.8. The guide field deviates the background H + 
in the outflow region: the H + outflow motion is tilted by the guide field and asymmetric 
with respect to the X point. The + acceleration is still localized along the separatrices, 
and normal to them. The peak velocities for H + and + with and without guide field are 
very close in value. 
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Figure 13 presents a comparison between the H + and + densities, normalized to no, 
at time Vl cH t = 21.8. Both H + and + ions tend to concentrate in the same regions in 
proximity of the separatrices. These enhanced density regions are represented in red in 
Figure 13; the depleted density regions are called cavities in the literature [Kleva et al, 
1995; Lapenta et al, 2010] and are visible in blue in the same plots. The formation of 
depleted and enhanced density areas along the separatrices is one of the signature of 
the guide field magnetic reconnection [Lapenta et al, 2010]. The maximum densities are 
almost identical for + and H + , while the minimum density for + is 3.6 times less than 
the one for H + . Figure 14 shows the peak densities at time VL cH t = 18.5: high density 
areas are not symmetric with respect to the X point, and are elongated thin structures in 
the case of H + , and are thick and shorter in the case of + . 

Figure 15 shows a different circulation pattern for the H + and + background species 
due the presence of a guide field. The guide field deflects the H + population in proximity 
of the separatrices in the inflow region and drives it to the X point. An inspection of 
panel b) of Figure 15 shows an + flow between separatrices: the + ions move from the 
cavities to enhanced density regions. 

4. Discussion and conclusions 

The simulations clearly show that the presence of an + population slightly decreases 
the reconnection rate and magnetic reconnection speed depends mainly on the lighter H + 
ion background species. This result is in agreement with previous Particle-in-Cell simula- 
tions [Hesse and Birn, 2004; Karimabadi et a/., 2010] with reduced mass ratios. However, 
simulations with physical mass ratio closely mimic real systems, providing realistic results 
that can be compared with observations. 
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It has been shown that magnetic reconnection with a H + current sheet occurs on the 

time scale of tens [Pritchett, 2001] in the presented simulations. Because the + ion 

gyro-period is 16 times slower the H + gyration, + particles rotate only once. The + 

population is essentially demagnetized over the time scale magnetic reconnection occurs 

in this simulation and its dynamics is determined by normal electric field accelerations 

in proximity of the reconnection fronts. This acceleration mechanism is already reported 

and studied in Refs. [Wygant et al, 2005; Drake et al, 2009b]. The average magnetic 

force exerted on + particles is typically four times lower than the one acting on H + ions, 

as shown in Figure 5. Differently from the + species, the background H + population is 

effectively magnetized out of the diffusion region and decouples from the magnetic field 

approximately at a du distance from the X point. The different dynamics of H + , + 

and electrons leads to a modified Hall current system, that produces a broader out-of- 

plane quadrupolar structure as reported in previous studies [Shay and Swisdak, 2004; 

Karimabadi et al, 2010]. The Hall magnetic field in presence of + shows a two-scale 

structure: an approximately 4d H thick region, where the B z (core Hall magnetic field) 

reaches the peak values, and an extended region, where B z (tail Hall magnetic field) 

decreases slowly. 

Because the H + and + background species undergo different acceleration mechanisms, 
that occur on different time scales and with different strengths, the two background species 
initially separate in the outflow region. The dipolarization front plasma is composed 
mainly of H + ions, while the + species follows the H + population in the outflow region. 

The reported simulations results lead to reasonable observational expectations in present 
and future missions. The tail Hall magnetic field is a possible signature of the background 
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+ presence during magnetic reconnection. Moreover, if H + ions are current carriers 

and a mixed background population of + and H + is present in the lobe plasma, the 

dipolarization fronts should be composed of a majority of H + plasma, and followed by a 

more extended + tail. 

Finally, the effect of a guide field on the + dynamics during magnetic reconnection has 

been studied. These simulations are of interest for reconnection in the magnetopause with 

+ ions [Pritchett, 2001; Lindstedt et al, 2010]. It has been found that the + presence 

does not change the reconnected flux and therefore the pace magnetic reconnection occurs. 

A comparison of simulation with and without guide field reveals that a guide field reduces 

the reconnected flux. The typical signatures of guide field magnetic reconnection, such as 

outflow asymmetry and the formation of region with enhanced and depleted regions along 

the separatrices, have been observed in presence of a background + also. In addition, 

an + flow, moving from the cavities to areas with enhanced density, has been identified 

for the first time. 

Two issues regarding to the validity of presented simulations should be finally addressed. 

First, it should be discussed if the grid resolution is fine enough to capture the electron 
layer physics and resolve the electron scales. Figure 16 shows the local grid resolution 
measured in terms of electron skin depth d e = c/u pe at time Q c nt = 17.5. The electron 
dissipation region is properly resolved (Ax/d e ~ 0.1). The local grid spacing Ax/d e 
is between 0.6 and 0.7 in the outflow area, where the species separation takes place. 
Thus, the electron scales are resolved by the grid in the regions of interest for this paper. 
However, the grid does not resolve the electron scales in the areas covered by exhaust 
plasma, where Ax/d e pa 2. 
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Second, the question if the + dynamics is allowed to properly couple to largest scales 

in a AOd H x 20d H box, should be addressed. In fact, the size of this domain is not 

large enough to accommodate the whole + dissipation layer, where the + ions would 

eventually couple to the H + species. This limitation is evident in Figure 4, where the 

quadrupole Hall magnetic field B z extends to the edge of the simulation box and vanishes 

because of the imposed boundary conditions. To answer this question, a simulation, whose 

domain is four times larger in each dimension, has been completed: a 160d H x 80d H box is 

discretized in a 4096 x 2048 grid and 6.5 billion particles are used. All the other parameters, 

physical mass ratios included, are the same of the simulation with smaller box. From the 

two panels of Figure 17, representing a pseudocolor plot of the Hall magnetic B z (panel a) 

and the B z profile along the line x = 68d H at time VL cH t = 17. 46 (panel b), it is clear that 

the larger simulation box includes the whole + diffusion region. The quadrupole Hall 

magnetic field, signature of the H + and + dissipation regions, extends over a limited 

area and vanishes in proximity of the current layer. The analysis of other quantities, that 

are not reported here, shows additional evidence that the 160<i# x 80d H simulation box 

includes the whole + diffusion region. 

All the features of magnetic reconnection in presence of + ions are recovered in the 

simulation with larger box. Magnetic reconnection starts later in time after approximately 

and proceeds at the same pace of the magnetic reconnection in the smaller box. 

Panel a) of Figure 18 shows a comparison of the reconnected fluxes of the two simulations. 

The bottom abscissa refers to the simulation with 40d H x 20d H box, while the top abscissa 

to the 160dn x 80d H simulation. A time lag of 1.5f2~^ is differentiating the two abscissae. 

The magnetic reconnection starts later in the larger box, the growth rate perfectly matches 
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in the two simulations with different box size and the reconnection rate does not change 

with the size of the simulation domain. 

Moreover, the separation of species occurs in the larger box simulation also. The panel 
b) of Figure 18 shows the peak densities for + and H + ions in the 160^ x 80d H sim- 
ulation box, confirming that the + and H + still separate in the outflow region, because 
of different acceleration mechanisms. In addition, the two-scale structure of the Hall 
magnetic field has been observed in the simulation with larger box also. Thus, after a 
comparison of the two simulations, it is clear that in the presented simulations there is 
no strong coupling of the + dynamics with the larger scales and that simulations with 
smaller box reproduce realistically the + behavior. 
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Figure 1. Evolution of the reconnected flux A\^, normalized to B cIh, and reconnection 
rate, normalized to the BqVa_h/c, for a simulation with only a background H + population 
with nbH = O.lno in black and for a simulation with mixed background population of H + 
and + ions with nbH = nbo = 0.05n Q in red. 
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Figure 2. Reconnection electric field E z , electric field X and Y components in units of 
{ujpHm H c) / e with superimposed magnetic field lines in blue at time VL cH t = 10.9. 
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Figure 3. Out-of-plane magnetic field component B z , in units of (u} p HmHc)/e, at time 
QcHt = 16.4 for a simulation with only a background H + population with nbH — O.Iuq in 
panel a) and for a simulation with mixed background population of H + and + ions with 
nbH = n bo — 0.05n in panel b). The quadrupolar structure reaches the same peak values 
(approximately 0.0009w p //m^c/e = 0.25£>o) in both cases, but it is broader in the case of 
the^smiulaW with mixed H+ Au l^ground 2 Bi^ula!t?6ns l . am D R 
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Figure 4. Hall magnetic field B z profile along the x = 8.0du line for the simulation 
with mixed + and H + (red line) and only H + (black line) background. The small peaks 
in proximity of the neutral line (y = 10dn) reveal the electron scales. In presence of an 
+ population, the Hall magnetic field B z profile (red line) presents a two-scale structure: 
sharp B z peaks (core Hall magnetic field) and slowly decreasing B z (tail Hall magnetic 
field). The B z peak values are approximately 0.0006u P HmHc/e = 0.16B . 
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Figure 5. v x B/c magnetic force magnitude for the H + and + background populations 
with superimposed magnetic field lines in blue at time Q c nt = 10.9. 
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Figure 6. In-plane average velocity for H + and + background species at time fi c #i = 
18.5. The arrow colors and lengths indicate the velocity magnitude, normalized to c. The 
maximum velocities 0.00461c and 0.001405c correspond to I.TJVah and 0.39Vah- 
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Figure 7. In-plane streamlines for the background H + and + populations in black 
and red with superimposed magnetic lines in blue at time VL c ut = 18.5. The black arrows 
show the flow direction. 
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Figure 8. Densities, normalized to n , for the background H + and + populations at 
time VL c ut = 18.5 with superimposed magnetic field lines in blue. 
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Figure 9. Peak densities regions for H + and + and electron background populations 
in black, red and green respectively at time VL c nt =18.5 with superimposed magnetic field 
lines in blue in panel a). In panel b) the density profiles for the three background species 
along the line y = L y /2 at time Q c nt = 18.5 with the same coloring scheme is shown. 
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Figure 10. Peak densities regions for H + and + background populations in black and 
red colors at three successive times. The H + and + populations progressively separate 
in the outflow region at the initial stage of magnetic reconnection. 
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Figure 11. Evolution of the reconnected flux A\l/, normalized to BqcIh for a simulation 
with mixed background population of H + and + ions with n&j? = n^o = 0.05rio, with a 
guide field in solid red line and without it in dashed red line. The reconnected flux with 
only background H + population and ti^h = 0.1n is presented in black. 
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Figure 12. In-plane average velocity of H + (panel a) and + (panel b) background 
species at time fl c nt = 21.8 in presence of a guide field. The magnitude of average 
velocity is normalized to c. The maximum velocities 0.004337c and 0.001582c correspond 
to 1.19V AH and 0.43 V AH - 
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Figure 13. Densities, normalized to no, for the H + and + background populations at 
time fl c nt = 21.8 with superimposed magnetic field lines in presence of a guide field. 
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Figure 14. Peak density regions for H + and + and electron background species in 
black, red and green with superimposed magnetic field lines in presence of a guide field 
at time Q c nt = 18.5. 
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Figure 15. In-plane streamlines for the background H + and + populations in black 
and red with superimposed magnetic lines in blue at time f2 c #£ = 21.8 in presence of a 
guide field. The black arrows show the flow direction. 
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Figure 17. Pseudocolor plot of the Hall magnetic field B z with superimposed black 
contour-lines to represent the non-zero areas of B z in the 160<i# x 80^ simulation box at 
time fl c nt = 17.46 in panel a). The line profile of the Hall magnetic field B z along the line 
x = 68dn is shown in panel b). The two plots reveal that the Hall magnetic field vanishes 
in a region close to the current sheet, and the whole + diffusion region is included in the 
lQOdn x 80cIh simulation box. 
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Figure 18. Comparison of the reconnected fluxes for simulations with AOdn x 20dn 
(proposed simulation box) and WOdn x 80dn (larger simulation box) sizes in panel a). 
Species separation still occurs in the 160d H x 80dn simulation box at time VL c ut = 17.5 
in panel b). 
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